Symmetry, complexity and multicritical point of the 
two-dimensional spin glass 

Jean-Marie Maillardf, Koji Nemoto$ and Hidetoshi Nishimori§ 

fLPTHE, Tour 24, 5 erne etage, case 7109, 2 Place Jussieu, 75251 Paris Cedex 05, 
France 

^Division of Physics, Hokkaido University, Sapporo 060-0810, Japan 
§Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro-ku, 
Tokyo 152-8551, Japan 

Abstract. We analyze models of spin glasses on the two-dimensional square lattice 
by exploiting symmetry arguments. The replicated partition functions of the Ising 
and related spin glasses are shown to have many remarkable symmetry properties 
as functions of the edge Boltzmann factors. It is shown that the applications of 
homogeneous and Hadamard inverses to the edge Boltzmann matrix indicate reduced 
complexities when the elements of the matrix satisfy certain conditions, suggesting 
that the system has special simplicities under such conditions. Using these duality and 
symmetry arguments we present a conjecture on the exact location of the multicritical 
point in the phase diagram. 



PACS numbers: 05.50.+q,75.10.Hk,75.50.Lk 
1. Introduction 

Properties of spin glasses are well understood in mean-field models, which show 
such anomalous behaviour as replica symmetry breaking, multivalley structure, and 
slow dynamics PJ E|. The difficult problem of whether or not these mean- field 
predictions apply to realistic finite-dimensional systems is still largely unsolved, and 
active investigations are carried out mainly using numerical methods [2]. Very few 
systematic analytical work exists, an exception being a symmetry argument using gauge 
invariance to derive the exact internal energy and several other exact /rigorous relations 

cum. 

In the present paper we develop another type of symmetry argument for models of 
spin glasses in two dimensions under the replica formalism. Our analyses reveal a variety 
of invariance properties of the replicated partition function under transformations of the 
edge Boltzmann factors, a notable example of which is the duality transformation. 

Also discussed are complexities of the edge Boltzmann matrix under inversions. It 
is well established that integrable systems such as the standard scalar Potts model have 
remarkably reduced complexities of the edge Boltzmann matrix when the parameters 
satisfy integrability conditions j^j- Analogous (but not the same) behaviour is observed 
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in our present problem, suggesting simplified (if not integrable) properties of the systems 
in restricted regions in the phase diagram. 

The results concerning symmetry properties are used to present a conjecture on the 
location of the multicritical point in the phase diagram. Although the argument is not a 
mathematically complete proof, we expect the prediction to be exact for several reasons 
including agreement with numerical results in many different models and satisfaction of 
inequalities. 

It should be remembered that most of the analyses are for real replicas, that is, 
the number of replicas is a positive integer. The quenched limit of vanishing number of 
replicas is discussed in relation only to limited cases including the conjectured location 
of the multicritical point. 

In the next section symmetries of the system are derived. The replicated partition 
function is demonstrated to be invariant under several different types of transformations 
of the edge Boltzmann factors. It is shown that a special subvariety exists which satisfies 
remarkable enhanced symmetries. Complexities under applications of matrix inverses 
are treated in section El Here also the same subvariety as above is seen to have reduced 
complexities, suggesting its special role. A conjecture on the multicritical point is 
presented in section |U and numerical evidence supporting the conjecture is discussed. 
The final section is devoted to summary and discussions. 

2. Symmetries of the partition function 

In this section we first derive the expression of the replicated partition function in terms 
of the edge Boltzmann factors. Then the symmetries and related properties of the 
partition function are discussed. The arguments are first developed for the ±J Ising 
model and then are generalized to a broader class of models (such as the Gaussian spin 
glass and random chiral Potts model) in the last part of this section. 

2.1. Replicated partition function 

Let us start by considering the Hamiltonian of the ± J Ising model 



where is ferromagnetic J(> 0) with probability p and antiferromagnetic —J with 
1 — p. The sum extends over nearest neighbours on the square lattice. The randomness 
in will be treated by the replica method. We will mainly consider the case of positive 
integer n, the number of replicas. Some results for the quenched limit n — > will be 
discussed in subsequent sections. The constraints of Ising spins and the ± J distribution 
of the interactions will be relaxed later in the present section. 




(1) 
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After the average over bond configurations, the system becomes spatially 
homogeneous^ , and the partition function 

[Z U ]a,v = Z n , (2) 

where the square brackets with subscript 'av' denote the configurational average, is 
specified uniquely by the edge Boltzmann matrix representing neighbouring interactions. 
For example, when n — 1 (the annealed model), the edge Boltzmann matrix is 



+ (l-p) 
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(3) 



W Wi 
Wi Wq 

where Wo = e K ,w\ = e~ K with K = J/ksT, and xo denotes the edge Boltzmann 
factor for parallel configuration of neighbouring spins whereas X\ is for antiparallel 
spins. Similarly, for n = 2 we have the edge Boltzmann matrix 
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This is a 4x4 matrix of hierarchical structure 

Xq X\ X\ X2 

X X Xo X 2 Xi 

X\ X2 Xo X\ 

X2 X\ X\ Xo 



A, B 1 
B 1 A, 



(5) 



with x = pWq + (1 — p)w\ 1 X\ = pw Wi + (1 — p)wiw and x 2 = pw\ + (1 — p)w^. The 
2x2 matrix B\ is obtained from A\ by replacing xq and x\ with x\ and X2, respectively: 
B\ = Ai(xo — > xi,xi — > X2). The element xo represents the two parallel spin pairs on 
neighbouring sites (such as ++ for replica 1 and ++ for replica 2), x\ is for one parallel 
and one antiparallel pairs (like ++, H — ), and X2 corresponds to two antiparallel pairs 
(H — , — h for example). 

The edge Boltzmann matrix of larger n can be derived by recursions. The general 
formula is 



P- 



Wo Wi 
Wl Wq 



+ (1-p) 



Wi w 
Wo w x 



A, 



A n -i B n _i 
B n -i A n _\ 



,(6) 



where B n _i = A„_x(xfc — > Xk+i',k = 0,1, ■ ■ ■ ,n — 1). The partition function © is a 
function of the matrix elements of A n , 

Z n {xo,X\, - ■ ■ ,x n ). (7) 

Here Xk denotes the Boltzmann factor of the spin configuration with k antiparallel spin 
pairs among n neighbouring pairs. 

\ except for boundary effects which are irrelevant to thermodynamic properties and will be ignored 
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In the symmetry arguments developed subsequently, we will often consider the case 
where the elements x k are independent of each other although they are originally related 
through the parameters of the ± J Ising model (K and p) by the relations 

Xq = p%Uq + (1 — p)w™ 

xi = pwQ^wi + (1 - p)w^ 1 w 

x 2 = pw^~ 2 wf + (1 - p)w T ^ 2 wl (8) 

x n = pw\ + (1 — p)wq . 
An advantage to consider generic points in the space spanned by the elements of the 
edge Boltzmann matrix 

S = {x = (x ,x 1 ,x 2 , ■ ■ -,x n )\x k E R (k = 0, 1, • ■ ■ ,n)} (9) 

is that we can discuss various well-known models such as the 2 n -state standard scalar 
Potts model which has Xi — x% — Xs — ■ ■ ■ — x n : All but one (xq) spin configurations 
have the same Boltzmann factors. The n-replicated ±J Ising model with the edge 
Boltzmann factor (jHJ) lies on a two-dimensional submanifold T of S: 

T ={xe S\x k = pw%- k wl + (1 - p)w^ k w k (k = 0, 1, • • • , n)} (10) 

2.2. Duality 

Some spin systems in two dimensions have invariance properties under duality 
transformations. The formulation of duality by Wu and Wang [6J is particularly 
useful for our problem having the edge Boltzmann matrix (JHJ) of hierarchical structure. 
According to these authors, the dual Boltzmann factors are derived simply by Fourier 
sums applied to each 2x2 block (corresponding to each replica) of the edge Boltzmann 
matrix (JEJ)§. The simplest case is the annealed model n = 1 of equation (jHJ): Its dual 
Boltzmann factors are the sum and difference of the original Boltzmann factors with 
appropriate normalization, 

V2x* = x + xi 
\f2x\ = Xq — x\. 

As an example, when the system is purely ferromagnetic p = 1 in equation (jH)), if we 
define the dual coupling K* by e~ 2K * = x\/xq (remembering e~ 2K = x\/xq = wi/wo), 
we have from equation the familiar duality relation of the ferromagnetic Ising model, 

e- 2K * = W °- Wl =tanhK. (12) 
w + Wi 

In the case of n = 2 with equations (j3J) and (jSj), it is necessary to generate combinations 
of sums and differences of appropriate matrix elements to obtain the dual Boltzmann 
factors, 

2x* = (x + xi) + (xi + x 2 ) = x + 2xi + x 2 

2x\ = (x - xi) + (xi - x 2 ) = x - x 2 (13) 

2 X 2 = (Xq - Xi) - (xi - X 2 ) = Xq - 2Xi + X 2 . 
§ Fourier sums for the two-component (Ising) case are just the sum and difference of two elements. 
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The formula for general n with the edge Boltzmann matrix (JBJ) is 



An = — 



1 

2~ 



' 1 1 


<8>n 


1 -1 _ 









" 1 1 


B n -i A n _i 




1-1 



(14) 



where the entries of the diagonal matrix A n are the dual Boltzmann factors 2 n l 2 x* m in a 
certain order. More explicitly, the dual Boltzmann factors, which are linear combinations 
of the original Boltzmann factors, can be written from equation ((Tl)) as 



2n/2 



JL .... 



k=0 



•^ki 



where the are the coefficients of the expansion of (1 — t) m ■ (1 + t) n 



(l-t) m -(l + t)' 



k=0 



that is, 



z=o 



-1^ 



n 



k 



m 
-I 



(15) 



(16) 



;i7) 



To understand equation (|TH|) we first note that m antiparallel pairs are chosen to generate 
x* m for dual spin pairs, which is reflected in the parameter m on the left-hand side of this 
equation. This is equivalent to choosing m of the second rows (1,-1) from the direct 
product of n 2x2 matrices in equation (|14j) ||. Then we choose k antiparallel pairs of 
original spins as in equation (JT7J) to obtain the coefficient of x^. 

On the square lattice the partition function remains invariant under the duality 
transformation of edge Boltzmann factors jHj: 

Zn{%Qi X X, %2, ' ' ' ; %n) = Z n [XQ, X*, X%, • ■ ■ , X*) (x E S) , (18) 

apart from a trivial factor 2 n and boundary effects, both of which are irrelevant to 
thermodynamic properties and will be ignored in this paper. Note that the symmetry 
under duality is valid for any values of the edge Boltzmann factors xq,xi,X2, ■ ■ ■ ,x n , 
which do not necessarily satisfy the relation (jHJ) of the ±J Ising model as indicated by 
the symbol x G S at the end of equation (|T8|) . 

It will be useful to write explicitly the general duality transformation (fTK|) applied 
to the specific case of the ± J Ising model for later use, namely, for the case x e T(c S). 
Again we first write the formula for the simple case of n = 2 written in equation (fE3|) 
so that the reader understands the structure: 



2x1 
2x\ 



p(w + Wif + (1 - p){wi + w ) 2 = (w + Wif 



Wo 



(19) 



= p(w + Wi)(w - wi) + (1 - p)(w 1 + w )(w 1 
-- {2p - 1) (w + w x ) (w - Wx) 

~- p(w - Wif + (1 - p){Wi - W ) 2 = (W - U7i) 2 . 

|| As seen in equation (|llf) . an antiparallel spin pair in the dual space corresponds to the difference 
(the row (1,-1)) of two original Boltzmann factors. 



2x% 
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The general formula is 

2 n/2 Z 2 *m = (w + W 1 ) n - 2m (w -W 1 ) 2m 

2 n/2 ^2m+i = (2p-l)(w + wi) n - 2m - 1 (u;o-wi) 2m+1 . 1 ' 

2. 3. Symmetries under sign changes of coupling constants 

In the present subsection the edge Boltzmann factors are regarded as functions of the 
parameters of the ±J Ising model following equation flSJ), x G T. The partition function 
Z n is invariant under the change of the sign of coupling constant K — > —K at all bonds 
if the system is on the square lattice. This change of the sign exchanges wq and w\, and 
consequently, according to equation (JSj), x k is exchanged with x n _ k - 

Z n (x Q ,x 1 ,x 2 , ■ ■ ■ ,x n ) = Z n (x n ,x n _i,x„_ 2 , ■ ■ • ,x ) (xeT). (21) 

Combination of two symmetries, duality ()18j) and change of sign of K (|21j) . leads 
to another symmetry. On the right-hand side of equation (|16j) the change of sign of t is 
equivalent to the change of the sign of D 1 ^ for odd k. The left-hand side of equation (fT7)|) 
suggests that the change of sign of t is also realized by exchange of m and n — m. Then 
we conclude in equation (|T5|) that the change of sign of Xk for odd k on the right-hand 
side 

D k m x k ^(-l) k D k m x k (22) 

should be performed simultaneously with the exchange of x* m and x* n _ m to keep this 
duality equation valid. Using duality (j!8|) we therefore find the following symmetry, 

Z n (%0, %2, —%3, ' ' ' 5 ( — l) n ^n) 

ry f * * 

— ^nK-^ni x n -\i ' ' ' i X \i X o) 

Zi n ^Xg, X-y, • • 4 , X n _i, X n ) 

= Z n (x ,xi, ■ ■ -,x n ) (xeT). (23) 

The second equality comes from equation ()21|) and the final relation is duality (fTHj) . 

When the number of replicas is even n = 2q, the partition function has additional 
symmetry which exchanges x k with x 2q _k for k odd only: 

Z 2 q(Xq, %1, %2i x ii " ' ■> x 2q-3i ^2g-2) %2q-l, %2q) 
= Z 2q (Xo,X2q-i,X 2 ,X2 q -3, X 4 , X 3 , X 2q _ 2 , X ± , X 2g ) (x E T) . (24) 

It is convenient to start the proof from the following expression of the edge Boltzmann 
factors generalizing equation (jHJ), 



xo = > Pie 



Xl =J2pie 2{q - 1)Kl 
i 

X2 =J2pie 2(q - 2)K ' 
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%2 g-2 = ^Pl e 



(25) 

-2 (g-2) if; 



%2q = y-pie 



I 



where the sum runs over 1—1,2 with pi = p,p 2 = 1 — p and K\ = K, K 2 = —K for 
the ± J Ising model. We would like the x k with even to remain invariant and the Xk 
with k odd to be changed into x 2q ^k- As one can see from equation (J25|) . this may be 
considered to be a transformation K —>■ —K for odd k only: 

x fc (AT) - x fc ((-l) fc -if). (26) 

The sign (— l) k actually corresponds to the sign of Ej • Ej for neighbouring sites i and 
j, where Ej = oj 1 ' ) af 2 ' ) • • • crf 2 ^ and Ej = cr^ crj 2 '* • • • crj 2 ^, because Ej • Ej is 1 for even 
number of antiparallel pairs and is —1 otherwise. The formula (J23j) can be written as 

x k = E^ ex p(^E^ a)f7 i a) V ( 2T ) 

J \ a=l / 

where k pairs among 2g pairs of neighbouring spins are antiparallel Yl a &i = 
2q — 2k. Changing K\ as 

K t -> (28) 

amounts to changing the edge Boltzmann factor as follows : 

x M -> J> exp ( if, • E, ■ E i £ a 4 w <xj Q) J . (29) 

/V a =l J 

In other words this is equivalent to performing for each replica a non-trivial Mattis 
transformation depending on all the other replicas: 

(a) (a) v ( a ) TT (£) /,n\ 

Note that, as far as dummy variables to be summed on are concerned, the r^'s are as 
good as the o~\ s : This is a one-to-one change of variables. Actually one goes back 
from the rf 1 ^ 's to the 's by performing the same transformation as the one defining 



the 7^ (a) 's : 



r («> _ a X«) = JJ T (ffl = H(^-af) = E^-E^K (31) 

We have used the identity E 2l? = +1. Therefore this Mattis transformation is an 
involution, and we have proved equation ()26|) which is equivalent to the symmetry ([24)1 . 

Let us point out another interesting symmetry of the partition function which can 
be derived from a combination of symmetries discussed so far. If we denote the duality 
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transformation (|15J1 symbolically as D and the exchange of Xk and x 2q _k for even n and 
odd k in equation ()24j) as M, the combination 

D 2 = D ■ M (32) 

also leaves the partition function invariant as long as the system satisfies the conditions 
for D and M to be the true symmetry such as the lattice structure (square lattice) and 
even n. This duality-like symmetry under the transformation D 2 will be used in the next 
subsection. It is worth noticing here that D and M commute for n even (D-M = M-D), 
and D, M and D 2 are all involutions, D 2 = M 2 = D\ = 1. 



2.4- Subvariety and duality 

The ±J Ising model with quenched randomness has remarkable properties along a line 
(curve) in the phase diagram defined by the relation 

exp(-2K) = — — (33) 
V 

known as the Nishimori line (NL) jHlHj- Similar interesting behaviour is observed on 
the same line also in the replicated system with n = 2 [2] . In this subsection we analyze 
symmetries of the replicated system with general integer n under the condition (}3*3*j) 
using the results of previous subsections. 

Let us consider the ±J Ising model on the square lattice whose parameters 
satisfy equation (}33|) . From equations (JEJ) and (}33|) . the latter being equivalent to 
p = wo/(wo + wi) and 1 — p — w\/{wq + wi), it is straightforward to see that the 
following relations hold: 

Xl = %n, %2 = Xn-1, ^3 = Xn-2, " " " , Xk = X n -k+l, • • • . (34) 

The condition ()33|) reduces the degree of freedom from two (p and K) to one. This 
means that we restrict ourselves to a one-dimensional curve in the (n + l)-dimensional 
space S. It will be useful to relax this constraint and consider the ([(n + l)/2] + 1)- 
dimensional submanifold of S specified only by equation (|3*4*|) . where [x] stands for the 
largest integer not exceeding x. We shall call this submanifold the subvariety J\f: 

N = {x e S\x k = x n -k+i (k = 1, 2, • • • , n)}. (35) 

The subvariety Af of course includes the one-dimensional curve NL specified by equation 
(IHSJ): NL C M. 

By the duality (fT3j) . M is transformed into the dual M* satisfying 

X\ = X 2 , X3 = £4, £5 = Xq, • • • , X2m-1 — %2m, ' ' ' • (36) 

To prove this fact, we first note that the condition ()34j) implies that the coefficient of 
Xk(— %n-k+i) on the right-hand side of equation (fl3j) is 

D k m (n)+D n m k +\n), (37) 
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where we have written the n-dependence of explicitly. Then the relation (J36)) follows 
if we can show 

DL-M + D^Wn) = D k 2m (n) + D n 2m k+ \n), (38) 

the left-hand side of which is the coefficient of Xk{= x n ~k+i) of 2 n ^ 2 x 2m _ 1 in equation 
([T3]) and the right-hand side is for Xk{= x n -fc+i) of 2 n ^ 2 x 2m . Equation (|55|) can be proved 
by induction with respect to n: The validity for small n(= 1,2,3) is checked trivially 
Let us assume that that equation (}3"Hj) is valid for n. Then it is not difficult to show 
that the same equation holds for n + l using the following recursion relation 

D k m {n+l) = D k m {n)+Dt\n), (39) 

which is derived by multiplying both sides of equation (jlfij) by 1 + 1 (which amounts to 
the change n — > n + 1). This ends the proof. 

Comparison of equations (pH} and suggests that the subvariety Af is in general 
not self-dual, Af ^ Af* . When n is even n = 2q, the partition function has an additional 
symmetry (|24j) . M in the notation of equation (J32|) . Then the combined transformation 
D2 = D ■ M keeps the subvariety invariant, Af = Af*: The application of D2 to Af is 
shown to yield the same relation as equation (J54*)l 

****** ** (a a\ 

"V» rf ry rf rf ry> ... r f* T' ... I /II I I 

x l — x 2<j! x 2 — x 2g-l! x 3 — x 2<j-2! 5 x m — x 2g-m+l) ■ 

The proof is outlined in |Appendix A We therefore conclude that the subvariety Af 



is globally self-dual when n is even, that is, any point in the submanifold defined by 
equation ()34|) is mapped by D2 to another point in the same submanifold (but is not 
fixed point-by-point in general). 



2.5. Inversions 

The edge Boltzmann matrix A n with generic values of the elements (x , x\, x 2 , ■ ■ • , x n ) 
has a remarkable property that its inverse matrix has the same structure; the inverse 
matrix has the same arrangement of elements as the original matrix. To show this 
we first note that the entries of the Boltzmann matrix are defined up to a common 
multiplicative factor, and therefore we can also discuss, instead of the matrix inversion 
A n — > A' 1 , a homogeneous matrix inversion transformation / (actually a homogeneous 
polynomial transformation): 

A n - A- 1 -det{A n ). (41) 

Let us also introduce the homogeneous transformation J corresponding to the inversion 
of elements of the dual matrix (Hadamard inverse), — ► l/yk (k — 0, 1, • • • , n) % 

m=n 

j ■ yk — > n Vm - ( 42 ) 

m=0,mj^k 

At first sight, performing the matrix inversion of the 2 n x 2 n Boltzmann matrix A n 
may seem to yield quite large calculations. However, since the duality transformation 

Transformation J is analogous to negating various coupling constants. 
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actually diagonalizes matrix A n as was mentioned in relation to equation (|14j) . it is 
straightforward to see that the matrix inversion just amounts to changing the dual 
variables x* m into their simple inverse : x* m — > ^/x* m - From this remark, it is 
straightforward to see that the matrix inverse of A n is a 2 n x 2 n matrix of the same 
form as the original A n (but of course with different elements : Xk — > x' k ). With 
obvious notations, and since there is no possible confusion, we will also denote by I this 
transformation on the x^s : 

J : x k — > x' k (43) 

Up to a common multiplicative factor one thus has (since D is an involution, that is, 
D 2 = 1) : 

D ■ I ■ D oc J, or I ex D ■ J ■ D. (44) 

Remark : One should recall that the two involutions / and J are actually 
(non-linear) symmetries of the phase diagrams of (anisotropic) spin edge lattice 
models [SJ |HJ E| • From equation (pf4*|) one sees that these "non-linear" symmetries of the 
phase diagrams are closely related to the (linear) duality symmetry. A duality symmetry 
exists when the edge Boltzmann matrix corresponds to cyclic matrices or semi-direct 
product of cyclic matrices |6 j . However, when a duality symmetry does not exist, like for 
instance the Ising model in a magnetic field, the two "non-linear" symmetries / and J 
still exist and can still be used to analyse the phase diagram [8J. 

2.6. Spin representation of the dual Boltzmann factor 

The elements of the dual edge Boltzmann factors of the replicated ± J Ising model (|2*U|) 
have an interesting symmetry. It is instructive to take the ratios of 
(which is equivalent to setting the energy level of all-parallel state to 0): 

(\ 2m— 1 
J = {2P ~ 1} tanh2m ~ 1 K (45) 

x* 2 Jxl = tanh 2m K. 

It is observed that the right-hand side is multiplied by tanh K each time the number of 
antiparallel spin pairs m increases. Another factor 2p—l appears alternately. Therefore 
the right-hand side of equation (J45|) can be expressed as a simple Boltzmann factor of 
the dual system, 

Aexp {K*(S W + 5 (2) + ■ • • + S in) ) + K* p S w S i2) ■ ■ ■ 5 (n) } , (46) 

where is the product of neighbouring dual spins in the oth replica (S^ = S^S^), 
and K* and K* are the dual couplings corresponding to thermal and randomness 
parameters: 

tanh K = e~ 2K * , 2p - 1 = e~ 2K i . (47) 

This expression (|4*U|) has a very interesting interpretation. The first part can be 
interpreted as being driven by thermal fluctuations since it has only the (dual) thermal 
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coupling K* in front of the spin variables. This first term is decoupled explicitly from the 
second part which is understood to be driven by quenched randomness as it is controlled 
by the (dual) coupling K* determined only by p. The first term causes ferromagnetic 
ordering of dual spin variables whereas the second term enhances spin-glass-like multi- 
replica ordering. 

The condition ()33|) is written in terms of the dual couplings as K* = K*. Therefore, 
on the NL, the above-mentioned two types of ordering tendency exactly balance. This 
reminds us of the result that the ferromagnetic ordering dominates above the NL whereas 
spin-glass order is larger below JU] • The balance K* = K* may also be called enhanced 
symmetry because the single-replica and the multi-replica terms have exactly the same 
coupling. 

It should be pointed out that a larger type of enhanced symmetry is realized by 
the 2 n -state standard scalar Potts model which has the edge Boltzmann factor 

Aexp {K^S^ + + ■ ■ ■ + S in) ) 

+ K^S^S® + S^S® + ... + fi-Cn-D^W) 

+ --- + K n S w S {2) ■■■S (n) } (48) 

with all couplings equal K\ = K 2 = K 3 = ■ ■ ■ = K n : It is straightforward to confirm 
that this Boltzmann factor has just two values, one for = = ■ ■ ■ = = 1 
and the other for all other configurations, thus representing the standard scalar Potts 
model. 

These two models, the ± J Ising model on the NL (with K* = K*) and the 2 n - 
state standard scalar Potts model (}4*%j) (with K x = K 2 = K 3 = • • ■ = K n ), coincide when 
n = 2 but not in general. We may therefore regard the n-replicated ± J Ising model on 
the NL as a system with similar but a little weaker symmetry than the standard scalar 
Potts model with 2 n states. Another point to notice is that the original Boltzmann 
factors (JHJ) of the ±J Ising model are not expressible in a simple form like equation (0HJ): 
If we try to do so by adjusting the coupling constants in equation (}48|) . each Ki will 
appear with a different value from the other coefficients. Only the dual system has the 
simple form of the edge Boltzmann factor (K 2 = K$ = ■ ■ ■ = K n -\ = 0). 

2. 1. Results derived from the dual Boltzmann factor 

The effective edge Boltzmann factor (|46J) shows that the replicated ±J Ising model in 
the dual representation has a relatively simple Hamiltonian with positive ferromagnetic 
couplings. This fact enables us to analyze the model using known results on 
ferromagnetic spin systems. One of our interests to be discussed later will be the 
structure of the phase diagram of the n-replicated system. The phase diagram drawn 
in terms of the original parameters p and T should look topologically the same as the 
one with axes 1/K* and K* because 1/K* is a monotone increasing function of p and 
K* is also monotonic in T (see figure Note that in the present paper we call the 
point marked in black dot in figure^ the multicritical point irrespective of the order of 
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transition. Of course one should take into account that ordered and disordered phases 
are exchanged between the original and dual systems. 







T, K* 


p 






SG 


/ F 





p = 1/2 p, 1/K * P =1 



Figure 1. Generic phase diagram of the n-replicated Ising model. The multicritical 
point (black dot) lies on the NL (shown in a thin dotted curve). 

We first check a few limiting cases. According to the dual Boltzmann factor (@EJ), 
when K* — > (p — > 1), the system decouples into n independent ferromagnetic Ising 
models as expected. If, on the other hand, K* — > (or the low-temperature limit T — > 
in the original variable), only the multi- replica coupling K* survives in equation (pfUj) . 
By redefining spins as Oi = S^S^ ■ ■ ■ S^ n \ the system reduces to a simple ferromagnetic 
Ising model. The critical point exists at K* = K^, where is the critical point of the 
ferromagnetic Ising model satisfying e~ 2Kc = a/2 — 1. Therefore the original system in 
the ground state has a critical point at p c = ^ = 0.7071 (derived from e~ 2K p = y/2 — 1) 
irrespective of n as long as n is a positive integer. Note that the quenched limit n — ► 
is expected to have a different critical probability p c near 0.89 (see the next section). 
This observation means that the two limits n —>■ and T — ► do not commute. 

The other limit K* — > oo corresponds to p — | in the original variable. According 
to equation (J46|) . the dual spins are then constrained as S^S^ ■ ■ ■ = 1. It follows 
that the nth spin variable can be expressed as = S^S® ■ ■ ■ S^" 1 ). Then the dual 
Boltzmann factor becomes 

Aexp {K*(S {1) + S {2) + ■■■ + S {n ^ + S {1) S {2) ■ ■ ■ S {n ~ 1} )} . (49) 

This is exactly the Boltzmann factor of the (n — l)-replicated system with K* = K* . 
We have therefore established the following relation for the partition functions of the n- 
and (n — l)-replicated ± J Ising models: 

Z n (K, K p = 0)<x Z n ^(K, K P = K), (50) 
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where the partition function is expressed by using the original couplings (K p defined by 
e -2K* _ tanh-fTp is the dual of K*). The trivial constants in front of the partition 
function are omitted. The identity (JHUj) proves that the n-replicated system with 
p — ~ (K p = 0) is equivalent to the (n — l)-replicated system on the NL, K = K p 
(equivalent to K* = K*). 

The effective edge Boltzmann factor ()46j) represents a system with ferromagnetic 
couplings only. Thus the Griffiths inequalities hold [TI]. In particular, first derivatives 
of an arbitrary correlation function are positive semi-definite: 

jL { S^sfS^ • • •> > 0, ^(St ] Sfsy • • •) > 0. (51) 

These inequalities imply that the boundaries between ferromagnetic and non- 
ferromagnetic phases are monotonic in K* and K* (and in the original variables T 
and p) as drawn schematically in figure ^ The reason is that, otherwise, the system will 
show reentrant behaviour and consequently the order parameters will be non-monotonic, 
violating (JSTj) . 

An interesting consequence of monotonic behaviour of the phase boundaries is that 
the transition temperature at p = ~ between the paramagnetic and spin glass phases (to 
be denoted as T c • (§)) is smaller than or equal to the temperature of the multicritical 
point Tc (MCP) which lies generically on the NL. Since we have already established 
in equation (JoDJ) that the n-replicated system on the line p — | is equivalent to the 
{n — l)-replicated system on the NL, it is concluded that 

Tj"-^(MCP) = T c (n) f~J < T c (n) (MCP). (52) 

In terms of p, this inequality reads, using K = K p (tanhK — 2p — 1) satisfied by the 
multicritical point to rewrite T by p, 

pf-^MCP) >p( n) (MCP). (53) 

The value of p at the multicritical point is a monotone decreasing function of n. We 
can then derive a lower bound for p^^MCP) of the quenched system n — > if we 
remember that the n = 1 annealed system is easily solved by reducing it to the regular 
ferromagnetic Ising model. The result is 



(MCP) > 1 + — - = 0.821797.... (54) 

Although this bound is not very tight since the expected value of the left-hand side is 
approximately 0.89 as discussed in the next section, it is nevertheless useful to have a 
mathematically rigorous bound (if we accept that the n — > limit causes no problems). 
We can also derive an inequality 

r w (I) < T^-D ( I) (55) 

from equation ()52j). 



2/ - c V2 
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Another non-trivial fact is the absence of reentrant behaviour as mentioned already. 
It is well established for n — 1 and n = 2 that the boundary between the ferromagnetic 
and non-ferromagnetic phases below the multicritical point is never reentrant as depicted 
in figure^in accordance with the present rigorous result j7j. It is a subtle matter whether 
or not this monotonicity still holds in the quenched limit n — - > 0: If this is the case, 
the phase boundary below the multicritical point should be a vertical straight line in 
accordance with the argument in [12] because the absence of ferromagnetic phase to 
the left of the multicritical point at any temperature is rigorously established in the 
quenched model |3]. However we avoid to make a definite statement here since it is 
not obvious that monotonicity derived from inequalities proved for positive integer n 
remains valid in the limit n — > 0. 



2.8. Relations between n and n + 1 

It has been established in equation (fBTij) that the n-replicated system with p = | is 
equivalent to (n — l)-replicated system on the NL on the square lattice. This result has 
already been pointed out by Georges et al p]. We here derive several further relations on 
physical quantities of n- and (n + l)-replicated systems. The arguments in the present 
subsection apply to the replicated ± J Ising model on an arbitrary lattice in an arbitrary 
dimension. No duality will be used. 

It is instructive to rederive the relation between the partition functions Z n and 
Z n+ \ without recourse to duality. The randomness average of the n-replicated partition 
function of the ± J Ising model is written as [3 H] 

Z "' g -^- (2co 3 hV ^ ey ' £r " Z(AT - (M) 

where N B is the number of bonds of the lattice, = ±1 is the sign of Jij(— ±</), 
and Z(K) is the partition function with fixed randomness. After gauge transformation 
and summation over gauge variables, the exponential in the above expression turns to 
a partition function with coupling K p [HE]: 

where N is the total number of spins. It readily follows from this equation that Z n (K, K) 
and Z n +i(K, 0) are essentially equal to each other: 

2^(2 cosh Kf B Z n (K,K) = 2 N+NB Z n+1 (K, 0). (58) 

This is the identity ()5U|) we already derived using duality for the square lattice. The 
present result is more general as we did not use the properties of a specific lattice. 
Similar identities hold for order parameters. The magnetization is defined by 

*»> = — Er ^,z ( A7. (59) 
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where (3 = 1/ksT. Spins on boundaries are fixed to up states to avoid trivial 
vanishing of the single-spin expectation value due to global inversion symmetry. Gauge 
transformation brings this equation into 

mJK, K p ) = Z " rVZ "' T — / \ >—_J. . 6 o 

This expression of magnetization is to be compared with the spin glass order parameter, 
which can be written using gauge transformation as 



n-2 



q n (K,K p ) - J2 T e KpJ:Tij Z(K)n 

E T (Es Sie-W) 2 Z{K p )Z{Ky-* 



J2 T Z(K p )Z(K)" 

Comparison of equations (JoTij) and (joTj) immediately leads to 

m n (K, K) = q n (K, K) = q n+1 {K, 0). (62) 

The first equality shows that there is no spin glass phase on the NL if we define the spin 
glass phase by m n = and q n > 0. The second equality of (Jo^j) indicates that the spin 
glass order parameter of the n-replicated system on the NL is exactly equal to that of 
the (n + l)-replicated system with p = \- 

The above identity (|62|) on the order parameters can be generalized to a relation 
between the distribution functions of order parameters [I]. We define and gauge- 
transform the distribution function of magnetization as 

P£\x;K,K p ) 



e K P Zn jZ (K) n 

E T (£s E. S(x 4E» Siaje-We-W*)) z{ky- x 



(63) 



E T Z(K p )Z(Kr 

where (3 p H(a) is the Hamiltonian with coupling K p and spin variables a. The 
distribution of spin glass order parameter, measuring the overlap of two replicas, is 
treated similarly. Using gauge transformation it is written as 

Pf\x-K,K V ) 



£ t e K P Y.n jZ (K) n 
E T (Es E. S(x 4E» SiaJe-WWM) Z(K P )Z(K)^ 



(64) 



J2 T Z(K p )Z(K)" 

It is easy to see from equations (|63|) and (J64j) that the following relation holds 

pW(x; K, K) = P( n \x; K, K) = P^ n+1 \x; K, 0). (65) 

Since the distribution function of a single- replica variable like magnetization P m (x) has 
only a trivial structure (with at most two delta peaks), it follows from the first equality 
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of (J65)) that the distribution of the spin glass order parameter P q (x) is also trivial on the 
NL for any n: There exists no complex structure such as replica-symmetry breaking. 
The second equality of (j65j) then proves that the same trivial structure holds for the 
spin-glass distribution function of {n + 1) -replicated system with p = \. Thus the spin 
glass phase of the (n + l)-replicated system with p = ~ has only a trivial structure. We 
should be careful when we consider the possibility to apply this result to the distribution 
function Pg n+1 ^ of the quenched system which, in the present context, corresponds to 
n + 1 or n ^ —1. The arguments in the present subsection are mathematically 
rigorous only for positive integer n. 



2.9. Generalization 

Many of the arguments so far have been on the replicated ±J Ising model although 
some of the symmetry relations are valid in more general systems with generic edge 
Boltzmann matrix with hierarchical structure. In the present subsection we apply the 
duality relations to the random Ising model with general coupling values as well as to 
the random Z q model which includes the random chiral Potts model. 

Let us first treat the n-replicated Ising model with a set of couplings K 1 , K 2 , K 3 , ■ ■ ■ 
with respective probabilities Pi,P2,P3, • • •■ The sign of Kj is arbitrary in this subsection: 
All the couplings can be ferromagnetic, for example. The hierarchical structure of the 
edge Boltzmann matrix is the same as the one already discussed for the ± J Ising model. 

The edge Boltzmann factors are 

3 
3 

X2 = Y,P3e (n - m > (66) 

3 

(67) 

e~ K >) 

e~ K i) 2 (68) 

2 n/2 < = EjPji^ ~e~ Kj ) n , 
a generalization of equations (fiT)|) and (|2()j) . It is sometimes convenient to take the 
continuum limit of randomness distribution. The original and dual Boltzmann factors 
are then 

Xk = /*,P(u) e *~»- (69, 



x n y Pj c j . 

3 

The dual are then 

2«/ 2 x* = J2j Pj ( eKj + e- K i) n -\e K i - 
2 n / 2 x* 2 = V . Pj \e K i + e -^) n - 2 (e^ - 
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2 n/2 x* m = I du P{u) (e^ + e - pu ) n ~ m {e Pu - e~ pu ) m : (70) 



A generalization of equation is expressed as the probability distribution of the 
form El 

P( u ) = e^Fiu 2 ), (71) 

where F(u 2 ) is an arbitrary even function of u and should be chosen to satisfy the 
normalization condition of P{u). The constraint (fTTj) applied to the edge Boltzmann 
factors (JoUj) and (fTUj) gives the symmetry properties 

X\ = x n , X2 = X n -i, £3 = x n -2, ■ ■ ■ (72) 

* =f= * * * * / r7n\ 

T* T* T* T* T* Y* ... 1 / -s 1 

x l — x 2' x 3 — x 4' x 5 — x 6' • \'°) 

The first equation (|72*j) can be checked directly from equations (jHTIjl and (|7T|l . The 
second (fTHj) has already been proved generically in equation When n is even, a 

further symmetry exists as in equation (|1D)) , implying that the subvariety A/" is globally 
self-dual, Af = Af*. 

We next turn to the random chiral Z q model for which again we use the Wu- 
Wang duality. The arguments in the previous subsections do not apply directly because 
we have used a binary (Ising) structure of basic variables. Each spin variable is 
now assumed to have q components to be symbolized by integers from to q — 1. 
The interaction will be written as V(k) when the difference in the neighbouring spin 
variables is k. The Ising model has q = 2, and the g-state standard scalar Potts model 
satisfies V{1) = V{2) = ■ ■ ■ = V(q — 1). The interaction is cyclic, V(k + q) = V(k). 
Chiral randomness is assumed to exist: The interaction energy for the difference in 
neighbouring spin variables being k is changed from the non-random value V(k) to 
V(k + I) with probability pi. 

The kth edge Boltzmann factor of the original lattice is 

x* = X>e BV < ,+ *>. (74) 
1=0 

The dual Boltzmann factors are obtained by Fourier transform of the edge Boltzmann 
matrix, a g-state generalization of equation (|14J) . and we only write the result here. See 
Appendix B for some details. 




q^x* =2Jwl2Je^l (75) 
1 \k=0 J 

q n ' 2 x\ =J2^{j2 eVik+l ^j (j2^e V(k+l ^j (76) 
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\ n—l 

( f^ x * q _ x =y,pi[ Y< eV{k+l) ) [ Y,^ l)k ^ v[k+l) ) (7* 



I \ k J \ k 



where u is the gth root of unity. These expressions, in particular those for x and x^, 
will be used in later sections. 



3. Complexity generated by inversions 

Since I and J are two (non-linear) symmetries of the phase diagram of the (anisotropic) 
model |H1 El , it is natural to combine these two involutions considering the infinite 
order (birational) transformation K = I ■ J. Keeping in mind that / is basically 
D ■ J ■ D as in equation f)44j) . we find that the iteration of K — I ■ J amounts to 
iterating k = D ■ J. 

The transformation K is generically an infinite order (birational) transformation 
which is a canonical symmetry of the parameter space of the model and an infinite 
discrete symmetry of the Yang-Baxter equations (star-triangle relations for spin-edge 
models) when the model is integrable This amounts to studying the iteration of 
transformation K and, in particular, the complexity of this transformation [T^j, namely 
the growth of the degrees ^3] of the (numerator or denominator of the) successive rational 
expressions encountered in the iteration 4- . Generically one gets an exponential growth 
of the calculations like X L with the number L of iterations. We will call from now on 
A the complexity of the transformation K or the complexity of the associated spin edge 
lattice model. In this "complexity" framework, an enhanced symmetry, occurring for 
some submanifold (in fact an algebraic subvariety) of the parameter space of the model 
(here the x^s), corresponds to a reduction of the complexity A to a smaller value : 
the integrable subcases pop out as algebraic subvarieties associated with a (at most) 
polynomial growth of the calculations [SHI] (A = 1). 

Let us first remark that the subvariety M is clearly (globally) invariant by J. The 
duality D permutes Af with its dual subvariety J\f*. The dual subvariety Af* (see 
equation (jUSl)) is also clearly (globally) invariant by J. Keeping in mind that / is 
basically D ■ J ■ D, one deduces immediately that Af is also (globally) invariant by /, 
and is, thus, is (globally) invariant by K = I ■ J (Af is (globally) invariant by k 2 where 
k — D-J). 

+ Or the purely numerical growth rate of the number of digits of the (numerators or denominators of 
the) successive rational numbers obtained when iterating an initial rational point [5]. Since the entries 
of the matrices grow exponentially, one has to use special representations of the integers allowing to 
manipulate large values (T^l 
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We have calculated the "complexities" A corresponding to the iteration of k = 
D ■ J, for a general 2 n x 2 n cyclic matrix, then for the general 2 n x 2 n matrix © 
for various number of replicas n and, finally, for Boltzmann factors restricted to the 
subvariety M of the 2 n x 2 n matrices © • The results are summarized in the following 
table : 



n 


n=2 


n=3 


n=4 


n=5 


Cyclic 


1 


5.8284 


13.9282 


29.9666 


Matrix © 


1 


2 


3.2143 


4.2361 


M 


1 


1.4142 


1.6180 


2.3344 



These values are all related to the solutions of polynomial equations with integer 
coefficients: They are the inverse of the smallest root of various polynomials with integer 
coefficients. For instance the complexities of a general 2 n x 2 n cyclic matrix, displayed in 
the first row of this table, correspond to the smallest root of polynomial 1— (2 ra — 2) t +t 2 
for various values of n. The complexities for generic 2 n x 2 n matrices of the form (JEJ), 
associated with n replicas, correspond respectively to 1 — t, 1 — 2 1, l — 3t — t 2 + t 3 , and 
1 — At — t 2 . Finally, the complexities for 2" x 2 n matrices of the form (jH} restricted to the 
subvariety Af, respectively correspond to 1 — t, 1 — 2 1 2 , 1 — t — t 2 , and 1 — 6 1 2 + 3 t 4 . 

We have also obtained the complexities A^ associated with the iteration of K — 
I ■ J, instead of A; = D ■ J . These A^-'s are (as it should, since K is equivalent 
to k 2 ) the square of the A's displayed in the previous table. For a general 2 n x 2 n 
non-cyclic matrix for which one does not have a duality anymore, one, surprisingly, 
finds, up to the accuracy of our calculations, the same complexities Xk as for a cyclic 
2 n x 2™ matrix ! The 2™ x 2™ matrices of the form (jBJ) are, of course, subcases of the 
general non- cyclic 2" x 2 n matrices. One thus finds quite a drastic complexity reduction 
((5.82842) 2 2 2 , (13.9282) 2 -> (3.214319) 2 , ...) when restricting 2 n x 2 n matrices 
to the form (JHJ) corresponding to effective Boltzmann matrix for n replicas. The replica 
analysis thus defines naturally a class of highly remarkable (non-random) lattice models, 
namely which are very interesting per se. We also see very clearly, that condition 
Af defines a highly singled-out subvariety of the previous remarkable models for which 
further drastic reductions of complexity occur (enhanced symmetry). 

It would be tempting to restrict models © to various (global or point-by-point) self- 
dual conditions (like Xq = xq, x* = x%, ...) in order to see similar drastic reductions 
of complexities, and possibly a polynomial growth of complexity on some integrable 
subvariety, or on some critical submanifold* (possibly given by intersections of the 
subvariety Af and various self-dual conditions) of the x , ■ ■ ■ x n parameter space S. 
Such an analysis, in fact, requires, on a square lattice for instance, to generalize our 
isotropic effective model (jBJ) to an anisotropic effective model with two different sets 
of horizontal and vertical Boltzmann matrices both of the form ©• The calculations 

* including the cases of first order transitions 
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become much more subtle and involved and, thus, will be detailed elsewhere. Let us 
briefly say that our preliminary "complexity analysis" results are not in disagreement 
with the numerical results displayed below. 

4. Conjecture on the location of the multicritical point 

We have discussed symmetry and complexity of replicated random spin systems and 
their generalizations. We now apply the results to the presentation of a conjecture on 
the exact location of the multicritical point in the phase diagram. 

4-1. General structure and conjecture 

In the homogeneous variables Xk, the duality relation we have derived for the ±J Ising 
model is (in the homogeneous partition function) 

Z n (xo, Xi, ■ ■ ■ , x n ) = Z u (xq, Xi, • ■ • , £*). (80) 

The question we try to answer in the present section is whether or not it is possible 
to identify the transition point with a fixed point of duality transformation for the 
replicated ± J Ising model and its generalizations (the Gaussian Ising model and random 
Z q model) as is the case for the non-random ferromagnetic Ising model and the standard 
scalar Potts model on the square lattice. Mathematically rigorous argument on this 
problem can be developed only in a few limited cases. We nevertheless take a step 
further to present a conjecture that the fixed-point condition x = Xq of only one of the 
n+ 1 independent variables has the possibility to lead to the exact transition point if we 
restrict ourselves to the one-dimensional subvariety with enhanced symmetry, the NL. 
Numerical data support our conjecture. 

Let us first consider the generic case of arbitrary values of xo,Xi,X2, ■ • ■ ,x n , not 
just on the NL, until otherwise stated. It is convenient to extract the factor xq, the 
edge Boltzmann factor for all-parallel spin state, in front of the partition function in 
equation (jHUJ): 

Z n (x , x 1 , x 2 , ■ ■ ■ , x n ) = (x ) NB Z n (ui, u 2 ,---, Un), (81) 

where Nb is the number of bonds on the lattice, and the reduced variables are defined 
by Uk = Xk/xo. This representation implies that we measure the energy relative to 
the all-parallel (perfectly ferromagnetic) state. The duality relation (}80|) is rewritten in 
terms of Z n as 

(x ) NB Zn{ui, u 2 , ■ ■ ■ , tin) = {x* ) NB Z n (ul,u* 2 , •••,<), (82) 
where u* k = x* k / Xq. 

When n = 1, this form of duality allows us to identify the fixed point with the 
transition point because the partition function has only one variable Z(ui), and u\ is a 
monotone decreasing function of u\. The transition point, if unique, should be at the 
fixed point U\ = u\. It then follows from equation (JH2*|) that Xq also satisfies the fixed 
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point condition x = Xq because the functional values of Z n on both sides are equal 
when the arguments U\ and u\ coincide. We can therefore use the fixed point condition 
xo = Xq in place of U\ = u\ to locate the transition point. If we apply this argument to 
the ±J Ising model (the annealed model), the condition xq = Xq defines a curve which 
is the exact phase boundary in the p-T phase diagram as depicted in figure 121 




Figure 2. Phase diagram of the ± J Ising model with n = 1. The phase boundary is 
given by xq = Xq and the NL is shown dotted. 

For n = 2, the situation is not very different. The condition xq = Xq is equivalent 
to Xi = x\ as well as to x^ = x* 2 as can be verified from equation (fH?j) . Thus we have 
Ui = u\ and U2 = u\ as soon as we impose the fixed point condition of the principal 
edge Boltzmann factor, xq = Xq. This fact is of course consistent with equation (|82j) 
because both the prefactor (xo) Nb and the functional value Z n (-ui,w 2 ) are the same on 
both sides when x — ^o- I n ^ ne case °f ^ ne Ising model, this fixed point condition 
Xq = Xq defines a curve in the phase diagram drawn in figure El (full curve above the 
multicritical point and dashed curve terminating at p = 0.5). As discussed by Georges 
et al j7], this curve coincides with the exact phase boundary above the multicritical 
point on the NL. Below the multicritical point the phase boundary splits into two parts, 
one separating the paramagnetic and spin glass phases and the other separating the 
ferromagnetic and spin glass phases. The duality transformation xo <-> Xq exchanges 
these two phase boundaries; the uniqueness assumption of the transition point is not 
satisfied in this case below the multicritical point. The NL has a high symmetry since 
it satisfies xi = xi as discussed in equation which is a globally self-dual subvariety 
because n is even, see equation pTJj) . This means that the NL is mapped to itself by the 
duality transformation and we are allowed to identity the transition point on the NL (the 
multicritical point) with the fixed point, x = Xq, which indeed gives the exact location 
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Figure 3. Phase diagram of the ± J Ising model with n = 2. The condition £0 = .Tq 
gives the correct phase boundary (full curve) above the multicritical point but not 
below (shown dashed). 

of the multicritical point. The lesson is that the condition x = x* does not always give 
the exact phase boundary but it may make sense to search for the multicritical point 
on the NL by investigating the intersection of the NL and the Xq = x* Q curve. 

Another straightforward case is the general 2 n -state standard scalar Potts model 
which has Ui = u 2 = u 3 = ■ ■ ■ = u n = u: The reduced partition function has then only 
a single variable Z n (u), and we can identify the duality fixed point u = u* with the 
unique transition point. It is worth pointing out that the condition xq = Xq together 
with equation (|82|) automatically means that the values of Z n on both sides are equal 
Z n {u) = Z n {u*) and consequently the arguments are equal u = u*. Therefore the fixed 
point condition of the principal edge Boltzmann factor x — ^0 sur h ces to identify the 
transition point (assuming it is unique) in the highly symmetric case of the standard 
scalar Potts model. Note that, when n — 2, the ±J Ising model on the NL coincides 
with the four-state standard scalar Potts model. 

Let us next move on to the case of general n. The subvariety Af defined as the set 
of points satisfying ui = u n ,u 2 = u n _i,u 3 = u n _ 2 , • • • is transformed into the dual Af* 
satisfying discussed in section EU In addition, the subvariety Af 

is globally self-dual A/" = A/"* (u\ = u* nl u\ = u* l _ l , ■ ■ ■) when n is even. The number of 
independent variables of Z n is [(n+ 1)/2] in Af, and it is not possible to straightforwardly 
apply the argument developed above for the simple case of a single variable. If we further 
restrict ourselves to the NL of the ± J Ising model, we consider the one- dimensional curve 
in Af representing the the condition which does not necessarily coincide with its 
dual (except in the simplest cases of n = 1 and 2). Therefore it is difficult to use duality 
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arguments to locate the transition point for generic n even if we restrict ourselves to 
the subvariety M (which satisfies a certain amount of symmetry) or its one-dimensional 
submanifold NL. 

We nevertheless proceed by trying an ansatz that the fixed point condition of the 
principal component of the edge Boltzmann matrix xq = Xq may lead to the exact 
location of the transition point if we restrict ourselves to the NL in the subvariety M . 
This ansatz is partly motivated by the above-mentioned example of the standard scalar 
Potts model, in which the condition x$ = Xq was sufficient to find the exact transition 
point. 

One may wish to try a different fixed point condition like x\ = x\ to investigate 
the possible location of the transition point. A physical reason to choose Xq = Xq, not 
the other ones, is that this edge Boltzmann factor xq is special in the sense that it is for 
the all-parallel spin configuration, and the expression of duality in equation ()82j) . which 
singles out Xq cLS £L special variable, means to measure the energy relative to the perfectly 
ferromagnetic state. This is natural when we discuss the properties of the system on 
the NL which starts from the perfectly ferromagnetic ground state (p = 1, T — 0) and 
terminates in the perfectly random high-temperature limit (p = |,T — > oo) and the 
multicritical point on this line separates the ferromagnetic and paramagnetic phases [3]. 
Another reason is that, if we try another condition, X ^ — X ^ for example, the result does 
not satisfy the rigorous inequality (JH2I)- This and other related points will be discussed 
in more detail later. 

The present ansatz xq = Xq to locate the multicritical point on the NL is equivalent 
to the one given in [T5] . 

The condition xq = Xq is written explicitly for the ±J Ising model as, using 
equations © and (1271 . 



This equation, derived from x — P ms * ne NL condition, gives the exact transition 
point (multicritical point) in the cases of n — 1 and 2 as already discussed. When 
n = 3, the curve given by equation ()83|) is drawn in figure 0] (marked 0) together with 
the curves coming from three other conditions and £3 = X3 marked 

1, 2 and 3, respectively. Apparently the intersections of the NL (dotted curve) and the 
four curves corresponding to be candidates of 

the transition point. To check these possibilities, let us note that equation coming 

(3) 

from Xq = Xq gives Tr = 1.65858 CIS db candidate of the multicritical point for n = 3. 
Similarly, = 1.75717 for n = 4, and = 1.82955 for n = 5 if we consistently 

use equation (|51j). When n = 1 and n = 2, the same equation leads to = 1.30841 

(2) 

and T c = 1.51865, both of which agree with the exact solution as they should. These 
conjectured values of from the condition x = x* Q plus the NL all satisfy the rigorous 



pe nK + (1 - p)e~ nK = 2~ n/2 (e K + e~ K ) n . 
This relation, in conjunction with the NL condition (jSHJ), yields 

e (n+l)K + e -(n+l)K = 2- n l 2 ( e K + e~ K ) H+1 . 



(83) 



(84) 




Figure 4. The curves representing x$ = Xq,X\ = x\,x% — x\ and X3 = x\ marked 0, 
1, 2, 3, respectively, for the n = 3 model. 



inequality (|52~j): 

T« < < < r c W < T c < 5 >. (85) 
If we instead assume that X\ = x\ together with the NL condition may give the 

(3) 

multicritical point, the same inequality is violated at n = 4 because, then, T c w = 1.56655 
and Tc^ = 1.51865. The same is true for the other Xk = x* k (k > 2). Thus xq = Xq is 
the only possibility consistent with the inequality if we are to choose the multicritical 
point among the intersections of Xk = x* k (k — 0, 1, • • • , n) and the NL. 

It should be of some interest to discuss the limit n — >• 00 because the problem can 
be solved exactly. The n-replicated partition function for arbitrary n is 

Z n = [Z(K) n ] av = [e- n ^ K \ v . (86) 

If we take the limit n — > 00 with N (the number of sites) kept large but finite, the 
number of terms in the above configurational average remains finite, and the value of 
this equation is dominated by the term with the smallest value of f{K). The smallest 
value of the free energy is expected to be given by the perfectly ferromagnetic bond 
configuration (and its gauge equivalents), and therefore we have 

Z n W e -PpNfo(K p ) e -nf3Nf (K)^ ^ 

where /o is the free energy of the non-random ferromagnetic Ising model. The first 
factor on the right-hand side is the probability weight of the non-random configuration 
and its gauge equivalents (proportional to Z(K p )\ see section |2~%|) . The system then 
has a unique critical point at the non-random critical point K = as well as at 
K p = (the latter being equivalent to p c = ^). The phase boundaries are two lines 
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representing these two critical parameters. On the other hand, the condition x = Xq 
gives in the limit n —>■ oo 

e - 2K = V2-l (88) 

which agrees with the above analysis giving K = Kj. Thus the conjecture that the 
multicritical point is given by the intersection of Xq = Xq and the NL is correct in this 
limit n — > oo. It is interesting to note that for n — > oo all the conditions Xk = x* k with 
k finite lead to the same conclusion. However, for finite n, the result depends on k, and 
the k = condition is the only one consistent with the inequality (J52j) . 

4-2. Numerical evidence for n = 3 

We have carried out an extensive Monte Carlo simulation of the n = 3 ± J Ising model 
on the NL on the square lattice to check the above-mentioned possible value of the 

(3) 

multicritical point T c = 1.65858. A preliminary scan of the energy showed that the 
transition is very likely to be of first order since clear hysteresis has been found. This 
is natural because the system has eight (= 2 3 ) degrees of freedom at each site and is 
not very far from the eight-state standard scalar Potts model (which has a first-order 
transition) because of the enhanced symmetry x\ = x^ on the NL (a little short of the 
high symmetry x\ = X2 = 23 of the standard scalar Potts model). 

We have therefore used the non-equilibrium relaxation method ( |17[ ITH] and 
references therein), which allows us to simulate very large systems in its initial relaxation 
stage to investigate equilibrium properties. Another reason to use the non-equilibrium 
relaxation method is that it allows us to identify the first-order transition point without 
evaluating the free energy ^7] . The results for the relaxation of magnetization are drawn 
in figure El for the temperature range from T = 1.650 (top curve) to T — 1.665 (bottom 
curve) with linear system size L = 8000, averaged over several samples, under mixed 
phase initialization appropriate for first-order transition jT2j. We have confirmed by 
comparison with smaller systems (and larger systems (L = 12000) in some limited cases) 
that this size L = 8000 is sufficiently large to investigate the properties of infinite-size 
systems at least to the time steps indicated in figure El (2000 steps). From the positive 
slope and upward curvature of magnetization after about 1000 Monte Carlo steps for 
T = 1.655, we conclude that this temperature is in the ferromagnetic phase due to 
the prescription of the non-equilibrium relation method J7j. Similarly the system is 
judged to lie in the paramagnetic phase at T = 1.665 from the negative slope and 
downward curvature. Our conclusion is Tc 3 ^ = 1.660(5), which is consistent with the 
above-mentioned conjecture 1.65858. The values from other fixed point conditions are 
completely ruled out, 1.56655 (from x\ = x{), 1.98207 (x% = x* 2 ) and 1.11466 (23 = x\). 

4-3. Quenched limit 

There have been a number of numerical studies of quenched systems (n — > 0) to find 
the location of multicritical point on the square lattice. If we take the limit n — > 
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0.50 



0.45 



0.40 



0.35 




Figure 5. Non-equilibrium relaxation of magnetization with mixed phase initialization 
for the n = 3 ±J Ising model on the NL on the square lattice. The temperatures are 
1.650, 1.655, 1.660, 1.663, 1.665 from top to bottom. 



in equation (jMj) . the formula rewritten in terms of p acquires a simple and appealing 
expression suggesting that exactly half of the entropy of bond distribution is exhausted 
at p c : 

In 2 



p c \np c - (1 - p c ) ln(l - p c 



9) 



Numerically this equation gives p c = 0.88997, which is compared very favourably with 
numerical results, among which 0.8894(9) is the most recent and extensive one [IB], as 
well as 0.8905(5) [H], 0.886(3) [20], 0.8906(2) [21] and 0.8907(2) [22]. 

Interestingly, the curve (|8H|) coming from xq = Xq is the only one with an intersection 
with the NL in the quenched limit. All the other curves Xk = x* k (k > 1) do not cross 
the NL as depicted in figure EJ 

If we further apply the condition xq = Xq to a model with continuous distribution 
of coupling as in equation ([69)1 . we find in the limit n —>■ 

J du P(u) ln(l + e~ Wu ) = (90) 

This expression already appeared in j2H]- Under the NL condition (fTT|) . this equation 
reads 

due Pu F(u 2 ) ln(l + e^ u ) = (91) 



which gives for the Gaussian model Jq/J = 1.02177. Numerical results are consistent 
with this value, 1.00(2) (Y Ozeki et al , private communication). 
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Figure 6. The curves representing fixed point conditions xq = Xq, x\ = x*, X2 = x\ 
and X3 = marked 0, 1, 2 and 3, respectively, in the quenched limit. The NL is 
shown dotted. 



P = jTTZ. 7 . ( 94 ) 



For the random chiral Potts model with 

7(0) = K, 7(1) = 7(2) = ■ ■ ■ = V(q-l) = (92) 
po = 1 - (q - l)p, pi = p 2 = ■ ■ -Pq-i = p (93) 

under the NL condition [24J 

1 

e K + q - 1 : 

the conjectured multicritical point coming from x = Xq is in the quenched limit, using 
equations (|71j) and (|7Hj) . 

- (1 - (q - l)p c ) ln(l - (q - l)p c ) - (q - l)p c lnp c = (95) 

Note that the model with equation (}9~3*j) reduces to the non-random case in the limit 
p —>■ in contrast to our previous convention of p — *■ 1; the reason is to avoid confusion 
when we compare our conjecture with the numerical investigation which suggests p c 
between 0.079 and 0.080 for q = 3 Equation (JHSJ yields p c = 0.079731 for q = 3. 
We therefore conclude that the condition xq = Xq has been confirmed to give the correct 
transition point on the NL within numerical accuracies for all available cases. 

5. Summary and Discussions 

We have derived symmetries of the models of the two-dimensional spin glass. In 
particular, invariances of the replicated partition function Z n under exchanges of edge 
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Boltzmann factors and duality have been proved. Also discussed is the invariance of the 
structure of the edge Boltzmann matrix by inversions. These properties hold for generic 
values of edge Boltzmann factors, not just for the replicated ± J Ising model. The dual 
Boltzmann factor of the ± J Ising model has an interesting expression in terms of dual 
spin variables. Griffiths inequalities apply to this expression, leading to monotonicity 
of phase boundaries, equivalent to monotonicity of the location of multicritical point as 
a function of n. 

Complexities under inversions have been investigated and a remarkable reduction of 
complexities (enhanced symmetry) has been observed in the subvariety TV. This result 
suggests that the behaviour of the system is simpler in this subvariety than at generic 
points even if the problem is not (generically) integrable because the exponential growth 
excludes integrability. 

Conjecture on the exact location of the multicritical point on the NL has been 
presented based on the duality and symmetry arguments. Reasons have been explained 
that the intersection of the curve xq = x^ and the NL is the only plausible candidate 
of the multicritical point among a set of similar conjectures (xk = x* k plus the NL). 
Numerical results are in very good agreement with this conjecture. Nevertheless, 
recalling the situation already encountered for the analysis of the phase diagram 
of the isotropic three-state chiral Potts model |2*71 |2*H] . where the critical manifold 
is numerically extremely close to the self-dual condition Xq = xo but is possibly 
mathematically different (in contrast with the situation encountered in the symmetric 
Ashkin- Teller model), one cannot discard the possibility that, restricted to the NL, the 
transition point could be numerically extremely close to the intersection with Xq = Xq 
but, actually, mathematically different. These points need to be further investigated 
carefully. 

It may be useful to consider the possibility that the condition xq = x* could be a 
very good approximation of the transition manifold away from the intersection with J\f 
but in the high temperature part of the phase boundary as in figure El at least for random 
ferromagnetic spin systems. For the g-state random standard ferromagnetic Potts model 
(the random version of the ferromagnetic scalar Potts model), the condition x = m 
the quenched limit is 



where the positive interaction K is assumed to be distributed with measure dfi(e K ). 
This is a Potts-generalization of equation (|HU|) . Investigation of the consequences of this 
equation is going on. 

Another interesting observation is that the formula gives K as a function of 
n which diverges as n approaches —1. It is not obvious at all that we can apply this 
formula to such a limit, but if we do so, then the divergence of K implies that the 
transition point on the NL is T c (n) -> as n — > — 1. If we use the relation between n 
and n + 1 in equation (|55|l . we find that the transition point of the n — > (quenched) 
system vanishes at p = \. Although we should be very careful to apply our results to 




(96) 
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such a limit of negative n, this last result is reasonable and interesting in its own right 
because it supports the usual consensus that the two-dimensional Ising spin glass does 
not have a finite-temperature phase transition when the distribution of bond randomness 
is symmetric. More efforts should be devoted to the investigation of this problem. 

Let us also comment here on the equivalence between the n-replicated 3<i random 
lattice gauge theory and its dual, the n-replicated Ising spin ferromagnet with edge 
Boltzmann factor described in section l2~o1 All the arguments sections 12 . 61 and f2 . 71 apply 
to such a case. In particular the bound on the critical probability like (|54jl results: 
^ 0) (MCP) > 0.9005. Numerically it is 0.97 |H1 jj. 

We have seen that the replica analysis naturally yields to consider a class of 
remarkable (non-random) lattice spin models, defined by equation (|5jl. which are 
highly structured and for which extensive symmetry analysis can be performed exactly. 
Recalling Domany's (N a , Np) terminology [2H], these models correspond to highly 
symmetric specific (A 2 , • • • , A 2 ) models and, more generally (N g ,---,N q ) models. 
These singled-out classes of lattice spin edge models provide a very powerful tool of 
analysis for many spin-glass problems and are also worth studying per se. 

One should note that many of the exact calculations, displayed in this paper, are 
still valid when the distribution of the coupling constants is not of the ± J type or some 
continuous distribution like equation (fTTj) but, for instance, a two-delta-peaks (Ji, J 2 ) 
distribution, or, even, a totally general distribution : We only need to get an effective 
Boltzmann matrix of the hierarchical form (JHJ). One should also underline that these 
models can straightforwardly be generalized to chiral Potts models (see | Appendix B ), 
and also to spin models without any Wu-Wang duality (like the Ising model with a 
magnetic field, see |Appendix B| ; the non-linear inversion relations I and J taking 
the place of the linear duality transformation D), providing room for many new exact 
results on spin-glass problems. 
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Appendix A. 

In this appendix we outline the proof of equation (J4"U|) assuming equation (|54*j) and the 
symmetry (f2~4"|) under the operation M. According to the duality relation ()15j) for even 
n = 2q, the expression of x* m , after applying the operation M, reads 

k=0 k=l 

jj This value is for the T = transition point, not for the multicritical point, but is likely to be close 
to the latter. 
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q g 

fe=0 k=l 

= D° rn x + ££x 2 + D^x 4 + ■■■ + D%x 2q 

+ D^x^! + D 3 m x 2q _ 3 + ■■■ + D%~ 1 x 1 . (A.l) 

We then impose the condition (jSlj) to find 

= D m x + (D l m + D 2 m )x 2 + (D 3 m + D 4 m )x, + - ■ ■ + (D 2 rr r 1 +D 2 ^)x 2q .(A.2) 

Thus, in order to show x* m = x 2q _ m+1 , it suffices to derive 

D 2k '\2q) + D 2k (2q) = Df q ^ +1 (2q) + D 2 2 k _ m+1 (2q) , (A.3) 

where we have written the n(= 2g)-dependence explicitly. This equation can be proved 
by induction with respect to q. 

The following relation will be useful for the proof: 

D k m+1 (n) + D k m - + \(n) = D k m (n) - D k m \n) (A.4) 

which is derived by replacing m in equation (|16|) with m + 1 (which amounts to 
multiplying both sides by (1 — t)/(l + t) ). 

It is easy to check explicitly using equation (fTTj) that the target relation (|A.3|) is 
valid for small q with any m and k. Let us then assume that this equation holds for 
q with any m and k and show that the same is true for q + 1. The left-hand side of 
equation (jA.3|) with 2q replaced by 2q + 2 can be reduced to an expression with 2q using 
the recursion relation (}39|) twice, 

Z^- 1 (2g + 2) + Z^(2g + 2) 
= D 2k {2q) + J D^- 1 (2g) + 2{D 2k -\2q) + ^(2?)) 

+ ^- 2 (2g) + Z^ 3 (2g). (A.5) 

Application of the other recursion relation ()A.4j) to each pair of terms on the right-hand 
side of the above equation yields 

^- 1 (2g + 2) + J D^(2g + 2) 
= D 2k _ x (2q) + Dlz\(2q) - D^(2q) - D^{2q). (A.6) 

This is our expression for the left-hand side of equation ()A.3j) with q —>■ q + 1. The 
right-hand side of equation (jA.3|) with 2q replaced by 2q + 2 can also be rewritten by 
the recursions ()39j) and (|A.4)1 to reach a similar expression 

D 2 2 k q +L m+1 (2q + 2) + Df q+2 _ m+l {2q + 2) 
= £ 2 2 J- m+2 (2g) + Df- l m+2 {2q) - D 2 2 k q Z 2 m+2 (2q) - D 2 2 k Z 3 m+2 (2q) . (A.7) 

This equation is equal to equation (jA.6|) by the starting assumption of induction, which 
completes the proof. 
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Appendix B. 



The results on the Ising model generalize straightforwardly to g-state models. Let 
us consider, for instance, the three-state chiral Potts model corresponding to a 3x3 
cyclic edge Boltzmann matrixff. For an arbitrary number n of replicas, the previous 
hierarchical scheme generalizes straightforwardly. If one denotes by A n the q n x q n 
effective edge Boltzmann matrices corresponding to [Z n ] av they can be obtained by the 
following recursion : 



A r . 



A 



n+l 





B n 


C 


C n 




B 




C n 


A 



(B.l) 



where 



A 1 



x 0,0 x 0,l x l,0 
x l,0 x 0,0 x 0,l 
x 0,l x l,0 x 0,0 



(B.2) 



and 



Br, 



A, 



rn.p 



X 



m,p+l)i 



Ar 



rn.p 



X 



m+l,p ) 



Note that the number of homogeneous parameters necessary to describe the pattern 
(IB.lj) grows quadratically with the number of replicas like (n + 2) • (n + l)/2 (this 
number is 3 for n = 1, 6 for n = 2, 10 for n = 3, 15 for n = 4, 21 for n = 5, 
...). The family of 3 n x 3 n models we have to study is thus a little more complicated 
(as well as the equivalent of the subvarieties M f!34j) ) . This is in contrast with the 
q = 2 Ising model without magnetic field, analysed in the text, for which the number 
of homogeneous parameters Xk grows linearly with n (this number is n + l). More 
generally, for a cyclic q x q matrix, the number of homogeneous parameters necessary 
to describe a pattern generalizing (jB.ljl will be the number of partition of n into q 



natural integers fcj : ki = k 2 + k 3 + ■ • • + k q = n with fcj = 0, 1, ■ • -n. The duality 



transformation on these models (jB.l|) is inherited from the duality transformation on 
the three-state chiral Potts model associated with cyclic 3x3 matrices. 

Let us consider a simple example of spin edge model with no duality transformation, 
namely the Ising model with a magnetic field. The effective Boltzmann matrices can be 
obtained similarly by a hierarchical scheme generalizing (jSJ) : 

A„ B„ 



where 



A 1 



A 



n+l 



x 0,0 x l,0 



B n C n 



(B.3) 



(B.4) 



^1,0 ^0,1 

ff That is the chiral Potts model studied by Baxter, Perk and Au-Yang 
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and 

Again the number of homogeneous parameters necessary to describe the pattern (|B.3j) 
grows quadratically with the number of replicas like {n + 2) • {n + l)/2. In this last 
case, we c?o no£ /iai>e a duality transformation D anymore : The inversion relations 
/ and J take the place of the linear duality transformation D. However, the inverse 
of matrices given by (jB.3|) are not matrices of the same form (jB.3|) . One thus needs 
to consider slightly more general matrices than (|B.3)1 in order to be able to use the 
non-linear symmetry /. Details will be given elsewhere. 
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